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ABSTRACT 



Context. Observations of dust polarized emission toward star forming regions trace the magnetic field component in the plane of the 
sky and provide constraints to theoretical models of cloud collapse. 

Aims. We compare high-angular resolution observations of the submillimeter polarized emission of the low-mass protostellar source 

NGC 1333 IRAS 4A with the predictions of three different models of collapse of magnetized molecular cloud cores. 

Methods. We compute the Stokes parameters for the dust emission for the three models. We then convolve the results with the 

instrumental response of the Submillimeter Array observation toward NGC 1333 IRAS 4A. Finally, we compare the synthetic maps 

with the data, varying the model parameters and orientation, and we assess the quality of the fit by a^ 2 analysis. 

Results. High-angular resolution observations of polarized dust emission can constraint the physical properties of protostars. In the 

case of NCC 1333 IRAS 4A, the best agreements with the data is obtained for models of collapse of clouds with mass-to-fiux ratio 

> 2 times the critical value, initial uniform magnetic field of strength ~ 0.5 mG, and age of the order of a few 10 4 yr since the onset 

of collapse. Magnetic dissipation, if present, is found to occur below the resolution level of the observations. Including a previously 

measured temperature profile of IRAS 4A leads to a more realistic morphology and intensity distribution. We also show that ALMA 

has the capability of distinguishing among the three different models adopted in this work. 

Conclusions. Our results are consistent with the standard theoretical scenario for the formation of low-mass stars, where clouds 
initially threaded by large-scale magnetic fields become unstable and collapse, trapping the field in the nascent protostar and the 
surrounding circumstellar disk. In the collapsing cloud, the dynamics is dominated by gravitational and magnetic forces. 
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1. Introduction 



individual objects: NGC1333 IRAS 4A - ISM: magnetic fields - stars: formation - magnetohydrodynamics 



Magnetic fields play an important role in the star formation 
process. Molecular clouds are expected to form dense cores 
through a combination of loss of magnetic and turbulent support. 
Eventually, a molecular cloud core overcomes magnetic support 
("supercritical" stage), and collapses gravitationally. The mag- 
netic field is then pinched and strengthened in the central re- 
gions of the core, and is expected to assume an hourglass shape 
(Fied ler et al. 1993; IGalli & Shut 11993115: INakamura & Ll 
120051) . 

Aspherical spinning dust particles tend to align their small 
axis parallel to the direction of the magnetic field. Thermal emis- 
sion from such elongated grains is thus partially linearly po- 
larized, with the polarization vector perpendicular to the mag- 
netic field. Consequently, the polarized emission is a good tracer 
of the magnetic field. To test the influence of magnetic fields 
we compare high-angular resolution observations of the polar- 
ized emission measured at submillimeter wavelengths toward 
the low-mass protostar NGC 1333 IRAS 4A with non-turbulent 
magnetohydrodynamic (MHD) models of molecular cloud cores 
threaded by an initial uniform magnetic field. This first step will 



help to (/) select the best models to describe the structure and 
evolution of low-mass cores, and, (ii) to better understand the 
importance of the physical processes involved in their formation 
and evolution. In a subsequent paper we will consider models of 
magnetized molecular cores formed in a turbulent environment. 

The low-mass protostar IRAS 4A is an ideal test site for 
models of magnetized cloud collapse and star formation. BIMA 
spectropolarimetric observations at 1.3 mm have detected and 
partially reso lved the polarizatio n in both the dust and CO (2- 
1) emission dGirart et all 1 19991) . showing hints of a hourglass 
morphology of the magnetic field. Recent polarimetric obser- 
vations with the SMA at 877 /mi with a resolution of 1!'3 
(390 AU) have shown that the magnetic field associated with 
the infalling envelope has a clearly "pinch ed" morphology o n 
a scale of a few hundreds AU (see Fig. 1 in lGirart et al.ll2006h . 
This morphology resembles the hourglass shape that is predicted 
by the standard theory of low-mass star formation in a collaps- 
ing core with a re gular magnetic field dominating the irregu- 
lar (tur b ulent) one dFiedler & Mouschoviasl 1 1993b IGalli & ShiiL 
Il993alfbt INakamura & Lil 120051). Ap plying the Chandrasekhar- 
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Fermi equation, iGirart et all d 1999b derived a magnetic field 
strength in the plane of the sky (POS) of fipos ~ 5 mG, corre- 
sponding to a mass-to-flux ratio of ~ 1 .7 times the critical value. 

iGoncalves et al.l (120081) compared the position angles in the 
plane of the sky of the polarization vectors determined by 
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IGirart etaL ( 2006]) with the inclination of ma gnetic field lines 
of ideal (iGalli & ShuL fl993allbh and non-ideal fshu et al.Ll2006h 



MHD collapse models. They found a good qualitative agree- 
ment for a source with < 1 M Q and a mass-to-fiux ratio of 
~ 2 times the critical value. The present work is a step for- 
w ard in the modelization and methodology with respect to that 
of lGoncalves et al" ] (120081) . 

This paper is organized as follows: in sections [2] and[3] we de- 
scribe the target source IRAS 4A and the selected MHD models, 
respectively. In section [4] we describe the synthetic map gener- 
ation. General results are detailed in section [5] In section [6] we 
present the MHD models prediction convolved with the SMA in- 
terferometer and compare them with IRAS 4A observations. In 
section Q we present the ALMA maps of the MHD model pre- 
diction. Finally, in section [8] we summarize the results and list 
the conclusions. 



2. NGC 1333 IRAS 4A 

The Perseus molecular cloud is an active low-mass star form- 
ing region, located at a distance ranging from 230 pc to 350 pc 



( Ridge et al 



Girart et al 



2006). For this work, we adopt the value of 300 pc 



2006) . In the souther n part of the reflection nebulae 



NGC 1333. Uennings et all (Il987l) were the first to identify the 
protostar NGC 1333 IRAS 4. ISandell et~ai1 (Il99ll) resolved the 
system into two different components, IRAS 4A and IRAS 4B, 
separated by ~31". They measured a luminosity of ~ 28 L Q (at 
350 pc, 11 Lq at 220 pc) equally shared between the two com- 
ponents. Subsequent interferometric observations have revealed 
further multiplicity: IRAS 4A is itself a binary system. The two 
com ponents IRAS 4A1 and IRAS 4A2 are separated by 540 AU 
(O. lLav et all[T995llLoonev et alJl2000t IGirart et al.ll2006h . 

This low-mass stellar system is in a very early stage of 
evolution. IRAS 4A and 4B are still embedded in an dense 
molecular and dusty envelope. ISandell et alj d 1 99 lb derived from 
submillimetric continuum single-dish observations a mass of 
~ 9 M Q . Subs equent interferom etric observations derived a 
mass of 1.2 M n (|Girart et alll2006l) for the compact component. 
iDi Francesco et aL ( 20011) detected infall motions from inverse 
p-Cygni profiles observed in H2CO (3i2-2n) and N2H + (1- 
0). Single-dish CO (3-2) observations reve aled a NE-SW well - 
collimated outflow arising from IRAS 4A (iBlake et all Il995l) . 
IChoil (120051) reports, through interferometric SiO (1-0) obser- 
vations, a highly collimated NE-SW outflow with a projected 
position angle of ~ 19°, and hints of a N-S outflow. The au- 
thor proposes that IRAS 4A2 is powering the main outflow while 
IRAS 4A1 would power the secondary one. 



2.1. New data 

For this work we gen erated new observa tional maps of IRAS 4A 
combining compact dGirart et all I2006I) . and sub-compact and 
extended (Ching & Lai, priv. comm.) configuration SMA data 
(see left panel of Fig. Q]), consisting of 8 h tracks in polariza- 
tion mode at 880 /mi. The data reduction was performed us- 
ing MIRIAD, while the imaging was done using GREG from 
the GILD AS package. To o btain the maps we used a weight- 
ing robust parameter of 0.5 dBriggsL Il995l) corresp onding to a 
beam of 1'.'24 x 1'.'12, slightly smaller than that of IGirart etafl 
d2006l) . Adding the sub-c ompact configurati on improved the 
map with respect to that of IGirart et alj (I2006I) : the sampling of 
the larger scales is better and allows a better characterization of 
the circumbinary envelope. In addition, the extended configu- 
ration data help in separating the emission arising from either 



the embedded compact sources or the circumbinary envelope. 
The combin ed continuum (Stok es I) map is in good agreement 
with that of IGirart et al.l d2006l) . atlhough the emission is more 
extended and has a sharper morphology. Furthermore, the po- 
larized intensity map covers a larger area, has a slightly higher 
intensity peak and a more defined morphology. 

The dust emission of IRAS 4A arises from the cold circumbi- 
nary envelope and from the warm circumstellar material around 
each protostar. Since the focus of this paper is on the morphol- 
ogy of the magnetic field in the circumbinary envelope, we have 
subtracted the contribution from the circumstellar component to 
the SMA visibility data. To do so, we first derived a map of the 
longest baselines (100-260 kX) corresponding to a beam of 0'.'70 
x (X'46. At these M,v-distances, the emission from the circumbi- 
nary envelope is resolved out, and the only contribution from the 
dust emissions arises from the circumstellar material (see cen- 
tral panel of Fig. Q] and Fig|6]l. Then, the clean components of 
this map were subtracted from the original visibilities. Finally, 
we obtained a new map of the circumbinary envelope using the 
resulting visibilities (see right-hand side panel of Fig. Q]). 

Table Q] shows the main parameters of the emission associ- 
ated only with the circumbinary envelope: peak position, r^st, 
RMS, S v , /P eak , FWHM, N Hl , «h 2 and mass. The integrated flux 
is 4.1 Jy, correspondi ng to a ma s s of .8 M©, both slightly 
smaller than those of IGirart et al.l d2006l) as we were able to 
isolate the envelope. The optical depth of the dust emission at 
880 yt/m imply that the observations trace very deep into the 
source. Therefore, neglecting scattering (see Sect. 14. U and as- 
suming an anisotropic radiation field, the polarized dust emis- 
sion is probably ori ginated in the alignment of dust grains to the 
magnetic field (see ^M^^ 

for a comprehensive review 
on this topic). The envelope hourglass morphology of the mag- 
netic field is more evident than in earlier data. A new feature 
is the double peak in polarized intensity. The map also shows 
a significant depolariz ation toward the so urce main axis, which 
was not as clear in the IGirart et al.l d2006l) map. This feature can 
be explained in ter ms of projections effect s intensified by beam 
smearing (see, e.g jGoncalves et aill2005l) . 



3. Theoretical models 

We compare the dust polarization map of IRAS 4A described in 
the previous section with the predictions of three mode l s of ma g- 
neti zed cloud collapse The models of IGalli & Shul (Il993allbh 
and lAllen et all (l2003allbh give the density profile and magnetic 
field distribution of an infalling envelope surrounding a low- 
mass star, the two mo dels differing ma inly in the choice of the 
initial conditions. The IShu et al.l d2006h model is similar to the 
previous two, but contains a parameter representing the spatial 
scale where the diffusive effects associated to an electric resis- 
tivity (assumed uniform) dominate the evolution of the mag- 
netic field. Therefore, our analysis is not able to test the theory 
of core formation from iniatially subcritical conditions by am- 
bipolar diffusion. This can only be accomplished by spatially re- 
solved Zeeman obser vations of molecular c loud cores and their 
surroundings (see e. g. ICrutcher et all 120 09). In a following pa- 
per (Frau et al., in preparation) we will analyze synthetic po- 
larization maps of protostellar cores extracted from numerical 
simulations of turbulent clouds. 
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Fig. 1. Left panel: IRAS 4A combining SMA sub-compact, compact and extended configurations with robust parameter set to 0.5. 
The stars mark the two compact sources (see central panel). Top-right corner shows the u,v range used. The synthesized beam is 
1'.'24 x 1'.'12. Contours show the dust emission at 880 pm in steps of 6<x, from 6<x to 96cr, where cr = 0.02 Jy beam -1 . The pixel 
map shows the polarized intensity (see scale on the right-hand side). Blue bars show the observed magnetic field direction. Central 
panel: IRAS 4A compact components. Legends are the same as in the left panel. The u,v range is restricted to the longest baselines. 
The synthesized beam is 0770 x 0746. Contours range from 6cr to 36<x. Right panel: IRAS 4A envelope emission (see section |2~TT i. 
Legends are the same as in the left panel. Contours range from 6cr to 48<x. 



Table 1. IRAS 4A: envelope continuum emission at 880 pm and 
derived parameters ( a ). 



o-(J2000) (") 


3h 29m 10.520s 


(5(J2000) ( b ) 


31° 13' 31". 12 


Tdust CO 


50 K 


RMSi 


20 mJy beanr 1 


RMSq 


2.5 mJy beam - ' 


RMSu 


2.5 mJy beam -1 


si 


4.1±0.4Jy 


"^peak 


1.03 ± 0.02 Jy beam" 1 


FWHM ( d ) 


1156 AU (3785) 


sT 


160 ± 16mJy 


peak 


38.2 ± 2.5 mJy beam -1 




14 arcsec 2 


T 


0.07 


Mt 2 (') 


1.2 x 10 24 cm" 2 


«H 2 ( f ) 


1.1 x 10 8 cm -3 


Mass ( f ) 


0.8 M Q 



( ) From a 2D Gaussi an fit to the source. 
n lGirartetaLl l l2006h . 

( d ) Diameter of the circle with the same area as the region of the source 
with intensity above half of the peak. 

( e ) Solid angle of the region with polarized intensity above 3-cr. 

( f ) Assuming k 2 sq qh z = 1-5 cm 2 g -1 and a gas-to-dust ratio of 100 
dGirart et alfl2006h . 



3.1. \Galli & Shil t1993afiR) 

This model follows the collapse of a singular isothermal sphere 
threaded by an initially uniform magnetic field. The cloud is as- 
sumed to be non rotating. This initial condition is a highly ide- 
alized representation of a non-equilibrium state. Inside the col- 
lapse region, bounded by an outward propagating slow magne- 
tosonic wave, the magnetic field dragged by the flow (even in 
the presence of ambipolar diffusion) deflects the infalling gas 



towards the midplane, forming a large pseudodisk. The initial 
state depends on two dimensional quantities, ro = 2a 2 / VGBo 
and to — 2a I "s/GBq, defining the characteristic spatial and tem- 
poral scale of the collapse. These depend on the sound speed 
a, the gravitational constant G and the initial (uniform) mag- 
netic strength Bo- F° r given ro and ?o, the time evolution depends 
on the non-dimensional parameter r = t/to, where t is the time 
elapsed since the onset of collapse. Fixing a = 0.35 km s , the 
model thus depends only on Bo and t. 

3.2. \Allen et all 12003M) 



This model is similar to that o f lGalli&Shul(fl993allbl) with some 
important differences: (/) being fully numerical, it overcomes 
the spati al and tempo r al limi tations of the semi-analytical ap- 
proach of lGalh & Shul dl993alfbl) : (ii) the initial state is a magne- 
tostatic unstabl e equilibrium con figuration (a "singular isother- 
mal toroid" see iLi & Shul Il996l) . already flattened in the direc- 
tion perpendicular to a magnetic field possessing a hourglass 
morphology from the start; {Hi) the cloud can rotate around 
an axis parallel to th e axis of the magnetic field. As in the 
iGalli & Shul (ll993albl) . magnetic field lines internal to a "sep- 
aratrix" are dragged into the accreting protostar. 

The initial configuration is specified by the sound speed, a 
(as in fhe lGalli & Shul[T993allb1 modeD. and the level of magnetic 
to thermal support, Hq, which represents the fractional overden- 
sity supported by the magnetic field above that supported by the 
thermal pressure, and the rotational speed, vo. The parameter Hq 
is related to the mass-to-flux ratio of the cloud. 

The flattening of the mass distribution (the "pseudodisk") 
and the magnetic field geometry are little affected by rotation. 
Conversely, the angular velocity of the infalling gas is strongly 
influenced by the magnetic braking associated to the strong field 
created by accretion, assuming ideal MHD. This effect has im- 
portant implications for the f ormation of rotati onally supported 
disks around young stars (see lGalli et all 120061) . 
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3.3. \Shu et all \200A) 

To overcome the difficulties associated to catastrophic magnetic 
braking and to the huge magnetic flux of the protostar, Shu et al. 

consider the consequences of non-ideal MHD effects dur- 
ing the accretion phase of low-mass star formation. In steady 
state, magnetic dissipation occurs inside a region of radius equal 
to the so-called "Ohm radius", rohm = ?7 2 /(2GM*), where 77 is 
the Ohmic resistivity (assumed uniform), G is the gravitational 
constant, and M* is the mass of the accreting protostar. Outside 
'"Ohm, the accreting gas is in free fall along radial field lines, that 
become straight and uniform inside rohm- The magnetic flux ac- 
creted by the central protostar is zero at all times. 

4. Synthetic map generation 

4.1. Assumptions 

To compare the intensity and the polarized intensity predicted by 
the models with the observed data, it is important to consider the 
effects of a temperature gradient, since the sub-mm emission is 
roughly proportional to the temperature. In this work, we have 
assumed both a uniform temperature profile (UTP) and a radial 
tempe rature profile (RTP) derived for IRAS 4A by iMaret et all 
(I2002h from water emission. Although the theoretical models 
considered here are computed assuming an isothermal equation 
of state, the IRAS 4A observed temperature gradient does not 
significantly affect the dynamics of collapse, because the kinetic 
energy due to thermal motions is more than one order of mag- 
nitude smaller than the kinetic energy of the inf ailing particles. 
For example, the temperature expected at 600 AU is 50 K, which 
leads to a thermal broadening of cr t henn~0.4 km s , whereas the 
infall velocity expected is Vff~1.7 km s _1 . 

We consider optically thin emission with no absorption or 
scatte ring effects, in agreement with the sub- mm emission prop- 
erties dHildebrandl 1 1 983b iNovak et al.L 1 19891) and with the opac- 
ity derived in IRAS 4A (see Table []]). We have assumed uni- 
form grain properties, represented by the polarizing efficiency 
parameter a which includes the a bsorption cross section and 
the alignment efficiency. Following iFiege & Pudritzl d2000l) . the 
maximum degree of polarization is 



Pmax — 



1 -a/6 



(1) 



We assumed p max = 15%, corresponding to a — 0.15. Despite 
the high value of a used, and the fact that the gr ain properties 
may change with density (IFiege & Pudritzl l2000h . we find that 
the absolute polarized intensities derived in the models match 
reasonably well the observed values in IRAS 4A. Lower values 
of a (e.g. a = 0.1) did not reproduce the data equally well. 

We performed the numerical integration using an equally 
spaced regular cubic grid and uniform step in the line-of-sight 
direction. 



4.2. Method 

We im proved the technique developed by iGoncalves et al.l 
d2008l) to compare theoretical models with observed data, in- 
cluding in the process the instrumental effects. In practice, we 
simulated all the steps of a regular observing run with both SMA 
and ALMA, generating synthetic maps with the same filtering 
and processing as the observed maps. With this technique we 
avoid any possible misinterpretation due to the effects of the in- 
strumental response and filtering, as well as the data modifica- 
tion because of the Fourier transform of the observed visibilities 



and the subsequent application of the dirty map cleaning algo- 
rithm. The process consists of a series of 5 consecutive steps for 
each realization: 



1. 



For any given model we generated three-dimensional (3D) 
data cubes of density and magnetic field components. The 
orientation in space of the 3D source models were defined 
by two viewing angles: the position angle (p of the projec- 
tion of the polar axis in the plane of the sky with respect to 
the north direction, and the inclination angle to of the polar 
axis with respect to the plane of the sky {to = 0° for edge-on 
view). Since the models used have axial symmetry, the op- 
tically thin emission assumption allowed us to explore only 
half of the inclination angle space (0° < a> < 90°). We re- 
stricted (f> to the range 0° < <p < 90° (the obs ervations fix 
the m agnetic axis of IRAS 4A at <p « 50°, see iGirart et all 
120061) . 

In the plane of the sky we simulated a square area with side 
length of 5 172 (~ 1 .5 x 10 4 AU at the distance of IRAS 4A). 
The map size was chosen to be about twice the SMA primary 
beam to better process the sidelobes in the final maps. In this 
plane we used a grid of 512 x 512 pixels with a pixel size 
of 071 (~ 30 AU), enough to oversample the smallest beam 
used in this work (~ 0.4"~ 120 AU for ALMA). With this 
choice we ensured, in the final convolved maps, the indepen- 
dence of points separated more than a beam distance due to 
beam convolution. In the line-of-sight direction we covered 
a length of 6 x 10 3 AU (equivalent to 20") sampled with 60 
cells. A larger integration length or a larger number of steps 
did not affect significantly the details of the final maps. 
Through a ray-tracing scheme, we integrated the emis- 
sion of the cells along the line-of-sight I generating 2D 
raw synthetic maps for the Stokes p arameter /, Q, and U . 
We followed a m ethod developed b y Lee & Draind d 19851). 
and el aborat ed bvlWardle & Konigll d 1 9901) . IFiege & Pudritzl 
d2000l) . and iPadoan et all d200ll) . One can calculate the 
Stokes Q and U intensities as Q — Cq and U = Cu, where C 
is a constant that includes all the terms assumed to be con- 
stant (polarization efficiency and polarization and absorption 
cross sections) that can be interpreted as a polarized intensity 
scale factor, q and u are the "reduced" Stokes parameters de- 
fined as 



-J 

u = j 



pB A {T&) cos 2ip cos y dl, 



pBx (7d) sin 2if/ cos y dl, 



(2) 
(3) 



where p is the density, B,i(7d) is the Planck function at the 
dust temperature 1// is the angle between the north direc- 
tion in the plane of the sky and the component of B in that 
plane, and y is the angle between the local magnetic field and 
the plane of the sky. Sto kes / is given by / = (C/a)(L - al.2) 
dFiege & Pudritzl 120001) where a is the maximum polarizing 
efficiency assumed to be 15%, the C/a factor can be inter- 
preted as a total intensity scale factor, while S and £2 are 
defined as 



-/ 



pB A (r d ) dl, 



// cos 2 y 1 \ 
pB A (T d ) (-^--1 dl, 



(4) 
(5) 
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and represent the emitted Stokes / intensity and the polariza- 
tion absorption losses respectively. As for Stokes Q and U, 
one can define the "reduced" Stokes / parameter as 

i = . (6) 

a 

4. Once the synthetic Stokes 2D maps were generated, the flux 
was rescaled so that the total flux in the region of the map 
with I v > 6cr (see section |2~TT i matched that of the synthetic 
map. A Gaussian noise was added to match the noise of the 
observations for each Stokes map. These maps were con- 
verted to visibilities using the same visibility coverage in the 
u, v plane as the real observations. In this step we mimicked 
the effects of the observation noise and the instrumental fil- 
tering. This allowed us to be sensitive to the same spatial 
scales and to similar emission levels. Then, the final model 
maps for the Stokes /, Q and U were obtained from the syn- 
thetic visibilities in the same way as the combined SMA 
maps presented in Fig.Q] 

5. The process described in points 1 to 4 is repeated a large 
number of times for both temperature treatments using dif- 
ferent values for (;) the position angle and the inclination 
angle a>, and (ii) the model parameters. 

For ALMA, we used the task simdata from the CASA pack- 
age to predict the expected maps for the models with the ALMA 
capabilities. At the source distance, the cells of the simulations 
had a typical length of ~ 30 AU. We chose the full ALMA con- 
figuration 09, which provides a synthesized beam of U.'l x 0'.'4 
(210 x 120 AU) thus ensuring that the final synthetic maps were 
not affected by resolution issues. We simulated a 2 hr run at 
345.8 GHz in polarization mode. As good weather is required 
for polarization measurements, we assumed 1 mm of precip- 
itable water vapor. The elevation of the source ranged between 
30° and 40°. 

4.3. Selection 

At this point, one has to select the "best" procedure for com- 
paring the synthetic maps with the observational data. This 
ca n be accomplish e d in s everal ways: (1) using the method 
of iGoncalves et al ] (120081) based on the minimization of the 
difference between the observed and predicted position angles 
of the polarization vectors (hereafter simply "angle difference 
method"); (2) performing a^ 2 analysis of the synthetic Stokes 
Q and U maps with respect to the observed maps. In the latter 
case (hereafter simply "x 2 method"), we positioned the peak of 
the synthetic map on the peak of the IRAS 4A envelope (see 
Fig. [TJ, and we compared the synthetic map with all the region 
of the observed map with intensity larger than 3cr (see Table [TJ. 
As we focus on the polarized emission, we define the best fitting 
models as those which minimized the sum^ 2 = x 2 q + Xij- Stokes 
/ was excluded since it shows considerable dependence on the 
assumed temperature profile. 

We illustrate the re sults for the two selection methods for 
the lGalli & Shul d!993allbl) models with RTP (see Figs. |2]and|3). 
Fig.|2]shows the standard deviation of the distribution of differ- 
ences in position angles as function of the viewing angles co and 
<p for the Galli & Shu (1993a,b) models described before. Only 
runs with average difference value of the position angles with re- 
spect to those of IRAS 4A lower than 15° are shown. Best- fitting 
models are characterized by the smallest values of the stan- 
dard deviation. The uniform distribution of results makes evident 



the low discrimination power of the angle difference method. 
Conversely, the x 2 method allows to perform a more significant 
selection of the best-fitting models. Fig.[3]shows a difference of 
more than one order of magnitude between bad- and well-fitting 
runs thus providing a higher discriminating power among all the 
runs. Therefore, in the rest of this paper, we adopt the^ 2 selec- 
tion method. 

Galli et al. (1993) - Radial T profile 
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Fig. 2. Pixel map: Standard deviation of the difference of the syn- 
the tic polarization vector s with respect to the observed ones for 
the lGalli & Shul(ll993albl) model with radial temperature profile. 
White pixels represent excluded models (see text). Grayscale is 
on the right-hand side of the panel. 
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Fi g. 3. Pixel map-. Sum of the^ 2 of the Stokes Q and U maps for 
the lGalli & Shul (Il993albh model with radial temperature profile 
with respect to the data. Contours: 10, 10 1/6 , 10 2/6 , . . . , 10 3 x 2 
levels. 



5. Results 

5.1. Orientation angles 

Stokes Q and U maps show a significant dependence on the ori- 
entation angles of the source, and can be used to constrain the 
viewing geometry. Fi g. [4] illustrates the diff erent emission pat- 
terns arising from an lAllen et al 1 (I2003albh source, with Hq 
0.12 s !. vo = 0, and t = 2 x 10 4 yr, after varying the orientation 
angles. Stokes Q and U maps are shown for all the combinations 
of position (4>) and inclination (oj) angles of 0°, 30°, 45°, 60°, 
and 90°. Conversely, Fig.|5]shows Stokes / maps which depend 
marginally on the position angle, and only for small inclination 
angles. 



4 = 30" 



4 = 45 6 



= 60" 




1H fr-l-M 



Fig. 4. Dependence of the synthetic Stokes Q and U maps on the orientation angles. The maps shown corresponds to the lAllen et al.l (l2003albl) model with H$ = 0.125, vo = 
and t — 2 x 10 4 yr. Thick contours: correspond to a single combination of position and inclination angle. Individual panels: Inside thick contours two panels are shown. Stokes Q 
and U map are shown in the left- and right-hand side panels, respectively. Common color scale for each Stokes map is shown below the first column. The angular scale is shown 
in the bottom right-hand side panel. Map contours: represent steps of 3-cr starting at 3-cr, where cr-2.5 mJy beam -1 . Columns: correspond to a position angle, 0, shown on the 
top. Rows: correspond to an inclination angle, a>, shown on the left-hand side of the figure. 
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Fig. 5. Same as Fig.|4]but for Stokes / maps. Each panel show the total emission (contours), the polarized emission (pixel map) and 
the magnetic field direction (segments). Map contours: for Stokes / represent steps of 6-cr starting at 6-cr, where <x=0.02 Jy beam -1 . 
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As can be seen in Fig.|H Stokes Q is specially sensitive to the 
position angle for small inclination angles, while Stokes U ap- 
pear to vary more with inclination angle at non-extreme position 
angles. From a practical point of view, one could identify with 
relatively high precision, by pure comparison with Figs. [4] and [5] 
both position angles of a source with a magnetic field with hour- 
glass morphology and an inclination angle smaller than ~ 60°. 
For inclination angles > 60°the expected magnetic field tend to 
be mostly radial and the Stokes Q and U maps show very sim- 
ilar morphologies independently of to. Therefore, polarization 
observations have shown to be a powerful tool to determine both 
position angles of a sources with respect to the LOS direction. 
Figures [4] and [5] could be used as templates for future observa- 
tions of the dust polarized emission toward star forming cores. 

5.2. Temperature profiles 

Figure|6]shows the differences in Stokes / u,v amplitude between 
both temperature treatments. The models show more realistic 
amplitudes using the RTP. The models with UTP show lower 
intensities than IRAS 4A in the > 25 kA range (< 4"), while the 
models with RTP fit the observations up to 15-20 kA (< 51'5). 
Note that this is remarkable for either temperature treatment as 
the intensity rescaling was done in the I v > 6<x which covered 
radii <3"(> 35 kA). Figs. [7J and [8] show the same information 
in the image domain for UTP and observed RTP, respectively, 
for cases with realistic model parameters and orientation angles. 
UTP maps show more extended emission (larger u,v amplitudes 
at short baselines) than IRAS 4A and a lower emission peak evi- 
dent from the residual maps. On the other hand, RTP maps show 
a slightly more extended morphology than IRAS 4A and real- 
istic intensity peak values. Note that RTP residual maps tend 
to show zero emission at short radii and slightly negative emis- 
sion at large ones due to the more extended sources predicted by 
the models. More realistic fluxes, intensity peaks and masses are 
derived from the RTP treatment. However, even using the ob- 
served temperatures, the radial intensity profile from models is 
steeper than that of IRAS 4A, which lead to smaller FWHM for 
model synthetic maps and, consequently, higher densities than 
those observed. 

The RTP treatment also predicted more realistic Stokes Q 
and U maps, shown in the middle left-hand side and middle 
right-hand side panels of Figs. [7j and [8] Although peak values 
are similar for both treatments, RTP map show roughly the same 
polarized flux over the same solid angle with similar morphology 
to IRAS 4A. On the other hand, UTP map showed roughly twice 
as much polarized flux over twice the solid angle of IRAS 4A 
(see Tables [2] and [3). An immediate consequence was the unre- 
alistic UTP vector map whereas the RTP one reasonably match 
that of IRAS 4A. 

Summarizing, RTP maps reproduced with higher fidelity the 
observed IRAS 4A emission, in the three Stokes parameters, bet- 
ter than UTP maps in the same conditions. Consequently, the 
physical parameters derived from RTP maps were more realis- 
tic. 

5.3. Visibility amplitudes 

Solid curves in Fig.|6]show the total u,v amplitude as a function 
of the u,v distance for the models using realistic model parame- 
ters and orientation angles (0=50° and a>=45°, see Sections [5. II 
and |6). The oscillations in the synthetic visibilities are due to 
the added noise. The intensity rescaling of synthetic data was 



performed in the image domain as it is the output of the sim- 
ulations. However, a good agreement in the u,v data is impor- 
tant given that it is the output from the telescope. Red and black 
dots in Fig. [6] show, respectively, the resulting u,v amplitudes 
before and after subtracting the two circumstellar compact com- 
ponents. We assumed that only emission from the envelope re- 
mained after subtraction. To test the goodness of the rescaling 
method we compared the observed u,v data of IRAS 4A with 
the synthetic u,v data derived from the models. A remarkable 
agreement (specially using the observed temperature profile, see 
Sect. 15.2b is achieved for u,v distances ranging from ~ 20 kA 
up to the maximum baseline with significant envelope emis- 
sion (~ 90 kA). None of the envelope models showed significant 
emission at u,v distances > 100 kA (equivalent to a radius of 1" 
or 300 AU) reinforcing the hypothesis that no envelope emission 
is detected from the IRAS 4A envelope in this u,v range. At u,v 
distances shorter than ~20kA the models show larger emission 
than the IRAS 4A envelope. At a distance of 300 pc this scale 
is equivalent to emission with a radius of > 1500 AU (> 5"), 
which is larger than the radius of IRAS 4A at a 3-cr level (~ 3", 
~ 900 AU). This excess has its origin in the fact that the mod- 
els predict a more extended source than the observed one, with 
typical radii of ~4"(~ 1200 AU) at a 3-cr level. 




20 40 60 80 '00 120 140 160 

u,v distance (kA) 



Fig. 6. Stokes I visibility amplitudes vs. u,v distance averaged 
in bins of 2kA. Top panel: Uniform temperature profile. Bottom 
panel: Radial temperature profile. Red and black dots: IRAS 4A 
full data and envelope data, respectively, with their statistical 
error bars. Vertical dashed line: u,v distance threshold used to 
derive the compact components map (see Fig. Q]). Starting at 
~90kA the emission seems to match with an unresolved source 
of ~1 Jy. Red, green, and blu e solid curve s : visibi lity amplitudes 
derived after c onvolving fhelGalli & Shul (1 1 993allbh . [Allen et al.l 
d2003alfbl) . and IShu et all d2006l) models, respectively, used to 
generate the maps of left-hand panels of Figs. |7jand[8] 
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Fig. 7. Comparison of selected models with data assuming a uniform temperature of the gas. The orienta tion angles are fixed to 
0=50° and w=45° for a bett er model comparison . Rows: IRAS 4A (first row); model lGalli & Shul (Il993allbh wit h B () _ 43 and 
r = 0.7 (second row); model I Allen et all (12003 albl) with H = 0.125, v = and f = 10 4 yr (third row); and model IShueTai] (120061) 
with rohm = 75 AU (fourth row). Columns: In each row, the panels show: intensity (first panel, contours), polarized intensity (first 
panel, pixel map) and magnetic field vectors (first panel, segments); map of Stokes Q (second panel, pixel map and contours); map 
of Stokes U (third panel, pixel and contours); residuals models-data for Stokes / (fourth panel, pixel map and contours). The color 
scale is shown on the top of each column. Contours: Contours for the Stokes I maps (left panels) depict emission levels from 6 cr 
up to the maximum value in steps of 6 cr, where cr = 0.02 Jy beam -1 . Coutours for the Stokes Q and U maps depict levels from 
the minimum up to the maximum in steps of 3 cr where cr = 2.5 mJy beam -1 . The solid red contour marks the zero emission level, 
solid white contours mark positive emission and blue dotted contours mark negative emission. Contours for the residual Stokes / 
follow the same rule of Stokes Q and U but with steps of 6 cr, where cr = 0.02 Jy beam -1 . The top right panel shows the beam and 
the angular and spatial scale. 

Table 2. Models with uniform temperature profile: 880 //m continuum emission and derived parameters ( a ). 





Unit 


Galli&Shu(1993a,b) 


Allen et al. (2003a.b) 


Shu et al. (2006) 


si 


Jy 


8.92 ± 0.25 


7.73 ± 0.25 


8.35 ± 0.25 


^Peak 


Jy beam 1 


0.77 ± 0.03 


0.79 ± 0.03 


0.98 ± 0.03 


FWHM ( b ) 


AU (") 


919 (3.06) 


1015 (3.39) 


892 (2.84) 




mJy 


390 ± 40 


301 ± 30 


328 ± 30 


.pol 
i Pcak 


mJy beam -1 


31 ±2 


34 ±2 


43 ±2 


n poT ( c ) 


arcsec 2 


28 


30 


23 


T 




0.27 


0.18 


0.30 


Mb ( d ) 


10 24 cm -2 


4.61 


3.14 


5.11 


»H 2 ( d ) 


10 8 cm" 3 


5.03 


3.10 


6.02 


Mass ( d ) 


M e 


1.92 


1.60 


1.82 



C) See Appendix A of lFrau e t al. (2010) for details. 

( b ) Diameter of the circle with the same area as the region of the source above the half peak intensity. 

( c ) Solid angle of the region with polarized intensity above 3-cr. 

( d ) Assuming k 2 so ghz=1-5 cm 2 g" 1 and a gas-to-dust ratio of 100 (jGirart et all l2006h . 
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Fig. 8. Same as Fig.Qfor the case with the lMaret et al.l d2002l) observational temperature gradient. 
Table 3. Models with radial temperature profile: 880 fim continuum emission and derived parameters ( a ). 





Unit 


Galli & Shu (1993a,b) 


Allen et al. (2003a.b) 


Shu et al. (2006) 


si 


Jy 


6.59 ± 0.25 


6.37 ± 0.25 


6.90 ± 0.25 


'peak 


Jy beam 1 


1.08 ±0.03 


1.02 ±0.03 


1.28 ±0.03 


FWHM ( b ) 


AU (") 


693 (2.31) 


744 (2.48) 


693 (2.31) 




mjy 


239 ± 24 


173 ± 17 


230 ± 23 


.pol 
*Peak 


mjy beam -1 


32 ±2 


30 ±2 


68 ±2 


Q. poi ( c ) 


arcsec 2 


16 


12 


12 


T 




0.37 


0.30 


0.39 


Nh, ( d ) 


10 24 cm -2 


6.28 


5.12 


6.64 




10 s cm" 3 


9.09 


6.91 


9.61 


Mass ( d ) 


M 


1.47 


1.39 


1.57 



(a) See Appendix A of Frau et al. (2010) for details. 

(b) Diameter of the circle with the same area as the region of the source above the half peak intensity. 

( c ) Solid angle of the region with polarized intensity above 3-cr. 

(d) Assuming k 2 so gh ; =1-5 cm 2 g" 1 and a gas-to-dust ratio of 100 dGirart et"aTl [2006), 
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6. SMA synthetic maps 

In general, synthetic intensity maps obtained from theoretical 
models tend to be less concentrated than the observed sources. 
The flux scaling based on the real data, combined with less com- 
pact synthetic sources, cause the synthetic maps to show more 
extended emission than the observed ones, and also a lower 
flux intensity peak. In the following subsections d6.ll 16.21 and 
\6.3\ we determine the best fitting parameters for each individual 
model (see Figs. [3] [9] andfTOli and perform a direct comparison 
of IRAS 4A with all the models (see Figs.|7]and[l and Tables|2] 
andE). 

6.1. \Galli& ShHt1993dB) 

We selected 5 values of the initial magnetic field (Bo = 0.86 mG, 
0.43 mG, 0.29 mG, 0.24 mG, and 0.17 mG, corresponding to 
f = 10 4 yr, 2 x 10 4 yr, 3 x 10 4 yr, 4 x 10 4 yr, and 5 x 10 4 yr) and 
3 values of the non-dimensional time r (r = 0.3, 0.5 and 0.7). 
The mass-to-flux ratio of the initial configuration is not spatially 
uniform as in the models of Allen et al. (2003a,b) described in 
Sect. 5.3. A spherical region centered on the origin and enclosing 
a mass M has a mass-to-flux ratio M/<p = nc 2 J{BoG 2 M). With 
the values of Bo listed above, and for a region enclosing a mass 
M = 1 M©, this corresponds to a mass-to-flux ratio, in units of 
the critical value, of 1.0, 2.0, 3.0, 3.6 and 5.1. 

For each choice of Bo and r we ran 70 cases corresponding 
to 10 different values of the position angle (p and 7 values of the 
inclination angle u>. We considered both an isothermal source 
and a radial temperature profile (see Section l4~Tb - For each of 
the 2100 maps generated, the model Stokes Q and U were com- 
pared with the observed values, and the sum of individual^ 2 was 
evaluated (see the radial temperature profile results in Fig. [3]). 
The results show that the best fit to the data is given by the mod- 
els with the highest values of the initial magnetic field Bq. In 
all cases, intermediate values of <p and o> are selected. Note also 
that for smaller values of Bo, the best fit is achieved for lower 
values of cj. This is due to a zooming effect: a larger Bo implies 
a smaller ro- For large Bo smaller angular distances mean larger 
radii, where the magnetic field configuration of the outermost 
parts of the model are naturally pinched in an edge-on view. On 
the other hand, for small Bo (large ro), the innermost region the 
magnetic field tend to be radial, and a larger inclination angle 
combined with the line-of-sight emission integration is needed 
to produce the pinched morphology. 

The fit is not sensitive to the value of the non-dimensional 
time r. Thus, the time elapsed since the onset of collapse is not 
well constrained by these models. For example, for models with 
Bo = 0.29 mG (third column in Fig. [3]) the time corresponding to 
the three values of r is 8.7 x 10 3 yr, 1 .5 x 10 4 yr, and 2.0 x 10 4 yr, 
whereas for models with Bo = 0.17 mG (fifth column), the time 
range corresponding to the three values of t is 1.5-3.5 x 10 4 yr. 

Fig. [7] and Fig. [8] show the predicted Stokes /, Q, U maps 
and the Stokes / residuals (second row) compared to the ob- 
served maps (first row) for this model, for the case with uni- 
form temperature and temperature gradient, respectively. Both 
realizations shown have Bo = 0.43 mG, r = 0.7 (corresponding 
to t - 1.4 x 10 4 yr after the onset of collapse), <p = 50° and 
u) = 45°. The derived physical parameters are shown in Tables|2] 
and [3] for the case with uniform temperature and temperature 
gradient, respectively. The radial temperature profile realizations 
show a better match to the observational maps for all the Stokes 
maps, as well as physical parameters closer to the observed to- 
ward IRAS 4A . The general morphology of the magnetic field 
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Fig. 9. Same as Fig. [3] for the I Allen et all (12003 allbh model with 
radial temperature profile. 



could be reproduced, as well as the depolarization toward the 
source axis and the double peak in polarized emission. 

6.2. \Allen et all 12003M) 

For the models of lAllen et all (12003 albh we selected 3 non- 
rotating models with different values of the parameter Ho defin- 
ing the mass-to-flux ratio of the initial state (Ho = 0.125, 0.25 
and 0.5, corresponding to a mass-to-flux ratio in units of the crit- 
ical value of 8.38, 4.51, and 2.66), and one rotating model, with 
Ho = 0.125 and uniform rotation velocity vo = 0.125 (in units of 
the sound speed). For each case, we considered 5 evolutionary 
times, t = 10 4 yr, 2 x 10 4 yr, 3 x 10 4 yr, 4 x 10 4 yr and 5 x 10 4 yr. 
As before, the position angles <p and oj were varied over a grid 
of 10 x 7 values, in both uniform and radial temperature profiles, 
generating a total of 2660 maps (due to numerical problems it 
was impossible to simulate the case Ho = 0.5 with t = 10 4 yr). 
The resulting x 2 of the comparison with the observed Stokes Q 
and U maps for the radial temperature profile cases is shown in 
Fig. [9] 

It is evident from the figure that better fits are obtained with 
lower values of the time elapsed since the onset of collapse, a few 
10 4 yr. The results are not very sensitive to the mass-to-flux ratio 
nor to the rotation of the initial configuration. There is a clear 
degeneracy between time and inclination angle with respect to 
the plane of the sky: a more concentrated field (a more pinched 
hourglass) can be obtained by letting the model evolve, or by a 
larger inclination of the magnetic field axis. For this reason, the 
region of minimum^ 2 moves towards lower values of o> at later 
times. 

Fig- E]and Fig.|8]show the predicted Stokes /, Q, U maps and 
the Stokes / residuals (third row) compared to the observed maps 
(first row) for this model, for the case with uniform temperature 
and temperature gradient, respectively. Both realizations shown 
have Hq = 0.125, v = 0, t = 10 4 yr, = 50° and u = 45°. The 
derived physical parameters are shown in Tables [2] and [3] for the 
case with unifor m temperature and tem perature gradient, respec- 
tively. As for the lGalli & Shul(ll993 a.b) models, radial tempera- 
ture profile realizations show a better match. For this model, the 



12 



Pau Frau et al.: Comparing star formation models with interferometric observations. I: MHD collapse models 



Shu et al. (2006) - Radial T profile 




15 30 45 60 75 90 
Position Angle (°) 



Fig. 10. Same as Fig.[3]for the lShu et al1(l2006l) model with ra- 
dial temperature profile. 



magnetic field morphology, source axis depolarization, and dou- 
ble peak in polarized e mission could be reproduced like with the 
iGalli & Shul (Il993albh models, although the polarized intensity 
derived was smaller in this case. 

6.3. \Shu et ali l2006\> 

To test this model, we varied the Ohm radius from 5 AU to 
150 AU. This parameter controls the size of the region where 
magnetic dissipation takes place and the magnetic field lines are 
almost straight. A total of 980 maps were generated, for 7 values 
of the Ohm radius, 10 x 7 values of the position and inclina- 
tion angles, and both isothermal and radial temperature profiles. 
Fig. [10] shows the x 1 of the comparison with the observed Q 
and U maps for the radial temperature profile case. The best-fit 
models ten d to have rph m in the ran ge 10-100 AU, as also found 
previously dGoncalves et alll2008l) . 

Fig. Q and Fig. [8] show the predicted Stokes /, Q, U maps 
and the Stokes / residuals (fourth row) compared to the observed 
maps (first row) for this model, for the case with uniform temper- 
ature and temperature gradient, respectively. Both realizations 
shown have rohm = 75 AU, (p = 50° and o> = 45°. The derived 
physical parameters are shown in Tables[2]and[3]for the case with 
uniform temperature and temperature gradient, respectively. As 
shown by the figure, this model fails to reproduce the double- 
peaked distribution of polarized intensity. This feature is associ- 
ated to strongly concentrated, almost radial magnetic field lines 
in the central region, at variance with the almost uniform field 
morphology produced by magnetic dissipation. Another charac- 
teristic of this model is the relatively high degree of polarization, 
not supported by the observations. 



7. ALMA synthetic maps 

The ALMA resolution used (0!'7x0r4 ~210x 120 AU 2 ) was cho- 
sen to have several projected cells of the simulations (~30 AU) 
inside each beam, allowing the best comparison possible with 
the models avoiding resolution effects. The ALMA sensitivity 
and u,v coverage is far much better than that of the SMA and, 
thus, allows (/) to map with a much higher fidelity the polarized 
emission, and (;;) to detect emission from a larger and fainter 
region. 

Figure [TTJ shows the ALMA maps for the models shown in 
Fig. [8] using the radial temperature profile. The level of detail of 
the convolved maps was very close to the original maps thus an 
almost perfect morphology reconstruction is possible for all the 



Stokes maps with the ALMA capabilities. Furthermore, the com- 
bination of sensitivity and resolution achievable with ALMA 
makes possible to extract usable information from the maps in 
a spatial range ~ 10 times larger than the resolution. We de- 
rived large and accurate polarization maps for all of the models. 
As Fig. [TTJ states, it will be possible with ALMA to reach reso- 
lution and detail levels which will allow to differentiate among 
different models, and to select those matching better the observa- 
tions. In order to make this result more evident we marked in the 
left-hand side panel of each model a red circle depi cting impor- 
tant distances related to the models. In the case o f IGalli & Shul 
(fl993allbl) and lAllen et all (12003115) red circles depict the loci of 
the isothermal collapse wave (r = c s t) which can be compared 
dir ectly to Figs. 2 an d 6 of the original papers, respectively. For 
the IShu el al. (2006) model, red circle depict the 10 rohm dis- 
tance, comparable with Fig. 4 of their paper. The high power of 
reproduction of ALMA encourages polarization observation to- 
ward all of the sources as a detailed modeling will be possible 
with the onset of this powerful instrument. 

8. Summary and conclusions 

The new data used in this work allowed to obtain a much bet- 
ter u,v coverage for the IRAS 4A region than in previous works, 
and, therefore, more reliable maps. In addition, the data added 
from different telescope configurations provided larger base- 
lines, to resolve the compact components, as well as shorter 
baselines, to better trace the extended envelope emission. This 
significant improvement allowed to separate the embedded com- 
pact sources from the diffuse envelope. A good u,v coverage was 
essential to perform reliable comparisons with models. To this 
goal, we developed a selection method, the so-called^ 2 method, 
with larger discriminat ing power than in previous studies (e.g. 
dGoncalves et alll2008l) l. 

The new data confirm that the source emission is optically 
thin, with no absorption or scattering, as expected for sub-mm 
emission. The opacity derived from the Stokes / map of the en- 
velope is negligible, implying that the maps trace very deep into 
the source. As the scattering appear to be negligible, the origin 
of the alignment of the dust grains is expected to be due to the 
magnetic field. 

Despite the complexity of the NGC 1333 star forming re- 
gion, MHD models of single star formation assuming quasi- 
static initial conditions and a uniform (or nearly uniform) mag- 
netic field show remarkable agreement with the observed char- 
acteristics of IRAS 4A, like intensity, polarized intensity, and 
polarization distribution. These facts suggests that the dust po- 
larization pattern resulting from the density and magnetic field 
distribution of non-turbulent models may apply even in less ide- 
alized initial conditions than normally assumed. Once the ori- 
entation angles are consistently determined, the comparison of 
the data with models of magnetized collapse indicate that a 
strong initial field is required (Bo larger than a few tenths of mG, 
see sect. 13.11 ), that the source is very young (a few 10 4 yr, see 
sects. [3~D and 13.21 ), and that the scale where magnetic dissipa- 
tion occurs is below the resolution of current observations (see 
sect. 13.31 ). However, with the current level of sensitivity and u,v 
coverage it is not possibl e to clearly discrimin ate a mong differ- 



ent collap se models. The lGalli & Shul dl993a|lbh andlAllen et al.1 
d2003alfbf) models fit better than the lShu et alj d2006l) model, but 



no selection can be done between the former two models. 

In general, the models predict sources with a more centrally 
peaked core and a larger less dense envelope with polarized 
emission less concentrated. A more refined dust grain treatment 
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Fig. 11. Same as Fig.[8]for ALMA configuration 09. Contours: Contours for the Stokes / maps (left panels) depict emission levels 
from 5 cr up to the maximum value in steps of 50 cr, where <r=0.2 mJy beam -1 . Positive contours for the Stokes Q and U maps 
depict levels fr om 5 cr up to t he max imum in steps of 10 cr, wher e <x=0.02 mJy beam 1 . Negative contours follow the same rule. Red 
circles: For the lGalli & Shul (1 1 993allbh andlAllen et al.l (12003 albh models dep ict the lo ci of the front of the isothermal collapse wave 
(r = c s t, see Fig. 2 of lGalli & ShuLll993al and Fig. 6 of lAUen et ail 12003 al) . For the IShu et al.l (l2006h model it marks the 10 r 0h m 
distance (see Fig. 4 of lShu et al.Ll2006l) . 



could help in a more realistic emission treatment. Current SMA 
observations of IRAS 4A clearly favor models with a temper- 
ature gradient. The total emission maps derived from models 
with a temperature gradient show the right peak value but larger 
fluxes, steeper profiles and more extended morphologies than 
IRAS 4A. On the other hand, they predict the right fluxes and 
morphologies of polarized emission but lower peaks and softer 
profiles. Conversely, polarization maps obtained with uniform 
temperature profiles are more extended than IRAS 4A maps (see 
Section 13^21 . Therefore, the inclusion of a realistic temperature 
profile cannot be ignored in modeling the sum-mm emission of 
low-mass protostars. 

An important result from the simulations of protostellar en- 
veloped threaded with an hourglass magnetic field is the pos- 
sibility of deriving the orientation angles of real sources from 
polarization measurements. Figs. |4] and show that up to incli- 
nation angles of ~ 60° it is possible to estimate both position 
and inclination angles from the Stokes Q and U maps. For larger 



inclinations the magnetic field tends to be radial and the Stokes 
maps do not show significant differences. Therefore, Figs.|4]and 
|5] can be used as templates for future observations of the dust 
polarized emission toward star forming cores. 

Another remarkable result is the good agreement in the u,v- 
plane of the observed and synthetic visibility amplitudes ob- 
tained assuming a temperature gradient, except for the short- 
est baselines (< 20 kA). The observed deficit of emission in 
IRAS 4A with respect to the models suggest a sharper density 
decrease at scales of ~ 1500 AU than predicted by theoretical 
models. 

The ALMA simulations have shown the capability of this 
new instrument to distinguish fine details even between mod- 
els of the same family. The methodology used in this work has 
proved to be a powerful tool to compare observations directly to 
theoretical models in a consistent way and avoiding instrumental 
effects. Future polarization measurements with the ALMA will 
provide real power to select the best models to describe the struc- 
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ture and evolution of low-mass cores and, consequently, to disen- Wardle, M., & Konigl, A. 1990, ApJ, 362, 120 
tangle the medium conditions and the physics ru ling the process . Zweibel, E. G. 1990, ApJ, 362, 545 
Upco ming radiative transfer codes like ARTIST dPadovani et al.L 
1201 lh will facilitate this kind of studies. ALMA data, together 
with powerful radiative transfer codes, may be used together 
with the technique developed in this work to extract as much 
information as possible from the data and to constrain the mod- 
els. 
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